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Abstract 

We propose a modified holographic dark energy (MHDE) model with the Hubble scale as the 
infrared (IR) cutoff. Introducing the infinite extra dimension(s) at very large distance scale, we 
consider the black hole mass in higher dimensions as the ultraviolet cutoff. Thus, we can probe 
the effects of the IR infinite extra dimension(s). As a concrete example, we consider the Dvali- 
Gabadadze-Porrati (DGP) model and its generalization. We find that the DGP model is dual to 
the MHDE model in five dimensions, and the ACDM model is dual to the MHDE model in six 
dimensions. Fitting the MHDE model to the observational data, we obtain that Q m0 = 0.269±£$$, 
^fco = O-OOS^q'q]^) and the number of the spatial dimensions is N = 4.78^^. The best fit value 
of N implies that there might exist two IR infinite extra dimensions. 
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I. INTRODUCTION 

From the cosmological observations such as Type la supernova (SN la) jl], cosmic mi- 
crowave background (CMB) j^j], and large scale structure (LSS) [3], there are strong evidences 
for dark energy (DE) which drives the accelerated expansion of the Universe. The most naive 
DE candidate is the cosmological constant (CC) introduced by Einstein. Although CC is 
consistent with the cosmological observations, there exists a fine-tuning problem [4]: CC 
((2.3 x 10~ 3 eV) 4 ) is extremely small comparing to the known energy scales such as the 



reduced Planck scale Mpi = 1/v8ttG (about 2.4 x 10 19 GeV) in general relativity and the 
electroweak scale (about 91 GeV) in particle physics, and we do not have any symmetry 
which can protect the tiny CC against quantum corrections (supersymmetry must be bro- 
ken above the electroweak scale). This fine-tuning problem is the greatest challenge in high 
energy physics. Also, there is a coincident problem Qj: why the DE and dark matter energy 
densities are comparable today since their evolutions are different as the Universe expands. 
Therefore, cosmologists.arrd particle physicists have proposed some other_DE modeler 



example, quintessence [5|, phantom |6j, fc-essence [7|, tachyon [8(], quintom j9j, hessence |l0| . 



Chaplygin gas [ill], Yang-Mills condensate [121 ]. etc. 



The DE problem might be a problem in quantum gravity [13] . However, we do not have 
a complete theory of quantum gravity right now. Fortunately, an important progress in 
the studies of the black hole theory and string theory is the proposal of the holographic 
principle [3], which may be considered as a fundamental principle of quantum gravity and 
then shed some light on the DE problem. Using the effective quantum field theory, Cohen, 
Kaplan and Nelson suggested that the quantum zero-point energy of a system with size L 
should not exceed the mass of a black hole with the same size, i.e., L 3 py < LMp X , where 
Pv is the quantum zero-point energy density Thus, the ultraviolet (UV) cutoff scale 
of a system is connected to its infrared (IR) cutoff scale. Applying this idea to the whole 
Universe, we can consider the vacuum energy as DE with density pde = Pv- Choosing the 
largest IR cutoff L which saturates the inequality, we obtain the holographic DE density 

p DE = Zc 2 Ml x L- 2 , (1) 

where c is an unknown constant due to the theoretical uncertainties and can only be de- 
termined by observations. Interestingly, taking L as the size of the current Universe which 
is the Hubble radius H~ l , one finds that the DE density is close to the observed value. 
However, Hsu [lj| pointed out that this yields a wrong equation of state for DE. 



To solve this problem, Li [17| chose the future event horizon of the Universe as the IR 
cutoff, and found that the model is a viable DE model. However, there exists an obvious 
draw back concerning causality: the event horizon is a global concept of space-time, and the 
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existence of event horizon depends on the future evolution of the Universe, i.e., the event 
horizon exists if and only if the Universe is accelerating. So, the original holographic DE 
model has presumed the existence of accelerated expansion. To avoid the causality problem, 
Cai proposed the agegraphic DE model in which the age of the Universe can be chosen as 
the IR cutoff 18(| . A new version of this model, which replaces the age of the Universe by 



the conformal age of the Universe, was suggested as well [19|. Moreover, Gao, Wu, Chen, 
and Shen proposed the holographic Ricci DE model where the average radius of the Ricci 
scalar curvature is chosen as the IR cutoff 20 ] . Thephenomenological consequences of these 
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models have been studied extensively 

In this Letter, to obtain the accelerating Universe and solve the causality problem, we 
propose a modified holographic DE (MHDE) model with the Hubble scale as the IR cutoff, 
and the UV and IR connection is modified by using the black hole mass in higher dimensions. 
As a concrete example, we consider the Dvali-Gabadadze-Porrati (DGP) models which have 
infinite extra dimension(s) 2a, |26j. The characteristic distance scale in such theories is 



called the crossover scale r c (for definitions, please see Section II.). At the distances that are 
smaller than r c , we observe four- dimensional gravity. While at the distances larger than r c , 
we observe higher-dimensional gravity. So at the IR energy scale which is smaller than r" 1 
(or say at the distance larger than r c ) , we need to consider the effect of extra dimensions. In 
particular, the mass of the higher-dimensional black hole is different from that of the four- 
dimensional one. Therefore, the UV cutoff scale in the MHDE model is different from that 
in the original holographic DE model. Moreover, our model provides a way to probe the IR 
infinite extra dimensions: the DGP model is dual to the MHDE model in five dimensions, 
and the ACDM model is dual to the MHDE model in six dimensions. Interestingly, with the 
Hubble scale as the IR cutoff scale, we can not only obtain the observed DE density with 
correct equation of state, but also avoid the causality problem. The best fit value indicates 
that there might exist two IR infinite extra dimensions. 

This Letter is organized as follows. We present our model in Section II. And we discuss 
the observational constraints in Section III. Our conclusion is in Section IV. 



II. THE MODEL 



First, let us briefly review the DGP model and its generalization 25|, |26(. The theory 
is a brane-world model which is embedded in a space-time with (asymptotically) flat in- 
finite n extra space dimensions. All the standard model (SM) particles are localized on 
the D3-branes, and the corresponding cut-off scale on the observable D3-branes can be the 
grand unification scale or higher. Also, the gravity spreads over the whole 4 + n space-time 
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dimensions. Thus, the action is 



S 



M 2+n 



J d*xd n yVGn + J d 4 x^7j (t 



C 



SM 



(2) 



where x^ and y l are respectively the coordinates for the Minkowski space-time and the extra 
space dimensions, G = \det (Gab)\ an d Gab is the metric for the whole space-time, g = 
|det((7 A1J/ )| and g^ u is the reduced metric on the D3-branes, M* is the high- dimensional Planck 
scale, T is the brane tension, and £sm is the SM Lagrangian 25j, [26| . The characteristic 
distance scale in such theories is called the crossover scale. For five-dimensional space-time 
with n — 1, the crossover scale is 



25] 



(3) 



And for six- or higher-dimensional space-time with zero tension branes, i.e., n > 2, the 
crossover scale is 
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Mi 



pi 



(4) 



Interestingly, at distance smaller than r c , i.e., r < r c , we observe the four- dimensional 
gravity, while at distance r > r c , we indeed observe high-dimensional gravity. Moreover, 



limeter gravity measurements [28] and the accelerator, astrophysical 
29j . the lower bound on the value of M* is M* > 10~ 3 eV. Also, 



from the data on sub- mi 
and cosmological data 

from the solar system measurements of Newton's law, one obtains the upper bound on 
M*: M* < 10 12 eV for five-dimensional space-time, and M* < 10 4 eV for six- and higher- 



dimensional space-time [271. l30l]. 



The idea of holography may be used to solve the CC problem 15|, [3l| . Cohen, Kaplan and 



Nelson proposed that for any state with energy E in the Hilbert space, the corresponding 
Schwarzschild radius R s ~ E is less than the IR cutoff L 15|. Under this assumption, a 



relationship between the UV cutoff p^ 4 and the IR cutoff is derived, i.e., L z p^ ~ L 15 ] 



Here the relationship M ~ R s between the mass and the horizon of the Schwarzschild black 
hole in four dimensions is used. In iV + 1 dimensions with iV = n + 3, the mass of the 
Schwarzschild black hole is [32J 



M 



(n - i)n N .i 



„N-2 
11 



(5) 



where 8nG 



N 



16nG N 

M}~ N . If we can see the effect of extra dimensions, then we have the relation 



and the DE density 



L 3 Pd 



Pd 



(N - 1)Q n _ LlN _ 2 
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(6) 
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where d is an unknown constant accounting for theoretical uncertainties. If we choose the 
Hubble horizon as the IR cutoff, then we get the MHDE density 

_ d(N-i)n N ^ N 

Pd ~ 16nG N H ■ (8) 
The mass of the Schwarzschild black hole is 

M = ^"^"^ Mf-Vg- 2 . (9) 

Thus, the high-dimensional black hole mass can be much higher than the Planck scale in our 
scenarios. However, to compare with particle physics theory, we consider the Lagrangian 
density. In our models, we have the dark energy density as follows 

Pd = ^zil^Mf- 1 ^" 5 . (10) 

Thus, for the model with one or two extra dimensions, the dark energy density is proportional 
to M 3 L _1 and M%, respectively. Thus, they are consistent with particle physics theory. 
However, for the models with three or more extra dimensions, the dark energy density is 
proportional to M^~ l L N ~ 5 with N > 6, thus, these models may have pretty large dark 
energy density, which is not consistent with particle physics. Therefore, we only consider 
one or two extra dimensions, i.e., the space-time dimensions are five and six, respectively. 
In particular, for five- dimensional and six- dimensional space-time, we obtain the observed 
DE density for M* ~ 10 7 eV and M* ~ 10~ 3 eV, respectively. Now the Friedmann equation 
becomes 

H 2 + — =—(p m + p r + p d ), (11) 

(1 O 

where the matter energy density p m = p m o(ao/a) 3 and the radiation energy density p r = 
p r0 (a /a) 4 . Substituting Eq. (1TUI) into Eq. (ITT]) , we get 

'i - a m0 - n k0 - n r0 ) — = n m0 (i + z f + n r0 (i + z f + n k0 (i + zf. (12) 



-"0/ V-no 



This is the same as the a dark energy model proposed by Dvali and Turner in Ref. 33] 
where a = 5 — N. The crossover scale in this model is 

w2/(JV-3) 

r (13) 



In [33j, the model is proposed as the modification to the Friedmann equation motivated 
by extra dimensions, and it was found that a < 1 to be consistent with astronomical 
observations. Here we use the holographic idea to derive the above Friedmann equation and 
take the model as a dark energy model. We also connect the holographic idea with the effect 
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of IR infinite extra dimensions. Furthermore, the parameter a has the physical meaning of 
the number of the spatial dimensions. If N = 4, i.e. , with one extra dimension, we recover 
the DGP model for the flat case. So the DGP model may be interpreted as the holographic 
dark energy model with the Hubble scale as the IR cutoff. When N = 5, i.e., with two 
extra dimensions, we recover the standard ACDM model. In other words, the cosmological 
constant is a special case of our model, and it may be interpreted as the effect of two extra 
dimensions. In our model, we choose the Hubble scale as the IR cutoff, so the problem of 
circular reasoning is avoided. In addition, the choice of the Hubble scale as the IR cutoff is 
more natural than other choices. 

Combining the Friedmann equation (fTTl) with the energy conservation, we get the equation 
of state parameter wa for the MHDE 

2a(l-n d ) a{Q m + 2Q r ) 
Wd + 3(2 -aQ d ) + 3(2 - aQ d ) ' 1 J 

So w d changes from —1 + 2a/ 3 during the radiation dominated era to —1 + a/2 during the 
matter dominated era. The deceleration parameter q = —a/(aH 2 ) is 

2n k + 3n m + 4n r 

?=-l+ -7Z • (15) 

I — aild 

Using the Friedmann equation Vt m + Q r + + tt d = 1 and q = 0, we obtain at the transition 
redshift z t 

(N - 3)Q k (z t ) + (N- 2)Q m (zt) + (N- l)Q r (z t ) = N - 3. (16) 

Since the transition from deceleration to acceleration happened very recently, we can ignore 
the contribution due to the radiation, and then we have 

(N - 3)fifco(l + zt) 2 + (N — 2)fi m0 (l + z t f = (N- 3)(H/H ) 2 . (17) 



III. OBSERVATIONAL CONSTRAINTS 



Now we use the observational data to fit the MHDE model. The parameters in the model 
are determined by minimizing x 2 = x 2 sn + yL p + xtmb + Xh- ^ OT ^ ne ^ a data, we use 
the Constitution compilation of 397 SN la 34|. The Constitution sample adds 185 CfA3 
SN la data to the Union sample 35j . The addition of Cf A3 sample increases the number of 
nearby SN la by a factor of roughly 2.6 — 2.9 and reduces the statistical uncertainties. The 
Union compilation has 57 nearby SN la and 250 high-z SN la. It includes the Supernova 
Legacy Survey 36j and the ESSENCE Survey 37|, [38(, the older observed SN la data, and 



the extended data set of distant SN la observed with the Hubble space telescope. To fit the 
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SN la data, we define 



v 2 

Asn 



307 

£ 

i=l 



[/iofes(^) ~ (J>(Zi)]' 

a? 



(18) 



where the extinction-corrected distance modulus n(z) = 51og 10 [di(^)/Mpc] + 25, /i G & s is the 
observed distance modulus, Oi is the total uncertainty in the SN la data, and the luminosity 
distance is 



d L (z) 



l + z 



:smn 



dz' 
eJz~' 



(19) 



where 



smn 



k\x 



em(y/\to k \x)/y/\QkU if < °» 
x, if fife = 0, 

k Sinh(-^/|fi fc |x)/^|n fc |, if fi fc > 0, 



(20) 



and the dimensionless Hubble parameter E(z) = H(z)/Hq. In particular, E(z) = [Q m o(l + 
zf + 1 - Ono] 1/2 for the flat ACDM model and E{z) = [a m0 (l + zf + (1 - fi m0 ) 2 /4] 1/2 + 
(1 — fi m o)/2 for the flat DGP model. We marginalized over the nuisance parameter Hq when 
evaluating x 2 sn - 

To use the baryon acoustic oscillation (BAO) measurement from the Sloan digital sky 
survey data, we define [39J 



r a (z d )/D v (z = 0.2) - 0.198 



0.0058 

where the effective distance is 



r s (z d )/D v (z = 0.35) - 0.1094 \ 
0.0033 ) 



Dv(z) 



dl(z) 



n l/3 



L(i + z) 2 #(*) 



(21) 



(22) 



The redshift ^ is fitted with the formulas 40j 

1291(fi m0 /i 2 ) - 251 



1 + 0.659(fi m0 /i 2 ) - 828 

6i = 0.313(fi m0 /i 2 )"°' 419 [l + 0.607(fi m0 /i 2 ) - 674 
and the comoving sound horizon is 

dz' 

r s {z) 



i + h(n b h 2 )% 

b 2 = 0.238(0^) 



(23) 

2^.223 ^24) 
(25) 



where the sound speed c s {z) = 1/y 3[1 + Rb/(1 + z)], the dimensionless Hubble constant 
h = H /100, and R b = 315000fi b /i 2 (T cmfe /2.7K)- 4 . 

To implement the Wilkinson microwave anisotropy probe 5 year (WMAP5) data, we 
need to add three fitting parameters R, l a and z#, so xlmb = Ax i Cov~ 1 (x i ,Xj)Axj, where 



c.Jz')E(z> 
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Xi = (R, l a , z*) denotes the three parameters for WMAP5 data, Axi = X{ — x° and 



Cov(xj, Xj) is the covariance matrix for the three parameters 4JJ- The acoustic scale I a is 



ird L {z.) 

La ~ (i + { ' 



where the redshift is given by 42J 



z* = 1048[1 + 0.00124(fi 6 /i 2 )-°' 738 ][l + g^Qmoh 2 ) 92 } = 1090.04 ± 0.93, (27) 

_ 0.0783(^/i 2 )~ a238 _ 0.560 

91 ~ 1 + 39.5(lVi 2 ) - 763 ' ^ 2 ~ 1 + 21.1(fi fe /i 2 ) 1 - 81 ' ( ' 

The shift parameter is 

R = 4= sinn f vW /" 4^) = L710 ± °' 019 - ( 29 ) 



z) at nine different red- 



via 

Simon, Verde, and Jimenez obtained the Hubble parameter H 
shifts from the differential ages of passively evolving galaxies 

thors in [44| obtained H(z = 0.24) = 79.69 ± 2.32, H(z = 0.34) = 83.8 ± 2.96, and 
H(z = 0.43) = 86.45 ± 3.27 by taking the BAO scale as a standard ruler in the radial 
direction. To use these 12 H(z) data, we define 



431 ] . Recently, the au- 



xl = j^ [H obs (z t )-H(z t )]^ (3Q) 
i=i 



where 07^ is the la uncertainty in the H(z) data. We also add the prior H = 74.2 ± 3.6 



km/s/Mpc given by Riess et al. 45j. The likelihood for the parameters in the model and 
the nuisance parameters Vt^h 2 and Hq (h) is computed using a Monte Carlo Markov chain 
(MCMC). The MCMC method randomly chooses values for the above parameters, evaluates 
X 2 and determines whether to accept or reject the set of parameters using the Metropolis- 
Hastings algorithm. The set of parameters that is accepted to the chain forms a new starting 
point for the next process, and the process is repeated for a sufficient number of steps until 
the required convergence is reached. Our MCMC code is based on the publicly available 



package COSMOMC |46j. 



By fitting the flat ACDM model to the combined SN la, BAO, WMAP5 and H{z) data, 
we find that \ 2 = 483.0 and f2 m0 = 0.272 ± 0.021. If we fit the observational data to the 
curved ACDM model, we find that x 2 = 482.9, fi m0 = 0.272^;°^ and tt k0 = 0.00 l+g;^. 
The contour plots and the probability distributions are shown in Figs. [T] and [2J By fitting 
the observational data to the flat MHDE model, we find that \ 2 = 482.5, tt m0 = 0.269l°:o25 
and N = 4.841q 40- The contour plots and the probability distributions are shown in Figs. [3] 
and|H Substituting the best fit values fl m o and N to Eq. fTTTI) . we get the transition redshift 
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A 2 



0.746. By fitting the observational data to the curved MHDE model, we find that 



482.2, a 



mO 



0.269 



+0.030 
-0.027' 



n k0 = 0.003±Sn" and N = 4.78 



-0.012 



+0.68 
-0.44- 



The contour plot and 



the probability distributions are shown in Figs. [5] and [61 By using the best fit values of Q m o, 
Jlfeo, and N, we find that the age of the Universe is to — 13.74 Gyr which is consistent with 
the result given by WMAP5 [41|. Substituting the best fit values of f2 m0 , flko, and N to Eq. 
(fIT]) . we get the transition redshift z t = 0.739. Since the best fit value of Qko is very small, 
the spatial geometry of the Universe is almost flat, so the transition redshift z t is almost the 
same for the curved and flat cases. In other words, we can neglect the contribution due to 
the curvature term. By using Eqs. (j!4p and (1151) with the best fit parameter values, we plot 
the evolutions of the fl m and fid, the equation of state parameter Wd, and the deceleration 
parameter q in Fig. 




0.24 0.25 0.26 



FIG. 1: The marginalized contours of fl m o and O^o in the curved ACDM model. 



IV. CONCLUSIONS 

We proposed the MHDE model with the Hubble scale as the IR cutoff, and our model 
fits the observational data as well as that of the ACDM model. The black hole mass in 
higher dimensions is used to modify the UV and IR relation, and then to derive our model. 
Since we used the Hubble scale as the IR cutoff, our model avoids the problem of circular 
reasoning. Furthermore, our model suggests a way probing IR infinite extra dimensions 
since the DGP model and the ACDM model are special cases of our model. In particular, 
the DGP model is dual to the MHDE model in five dimensions, and the ACDM model is 
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0.24 0.26 0.28 0.3 -0.01-0.005 0.005 0.01 

FIG. 2: The probability distributions of f2 m o and £lko m the curved ACDM model. 

dual to the MHDE in six dimensions. In other words, the effect of one extra dimension is 
manifested as the DGP model, i.e., the effect of five-dimensional gravity at large distance 
can be seen as DE in four dimensions. Also, the effect of six-dimensional gravity at large 
distance can be considered as the CC in four dimensions. Fitting the model to the combined 
SN la, BAO, WMAP5 and H{z) data, we find that Vt m0 = 0.269+°;°^, n k0 = 0.003±g;g^ 
and N = 4. 781° 44- By using the best fit values of O m o, Qko, and N, we find that the age of 
the Universe is to = 13.74 Gyr. The best fit value of N suggests that there may exist two 
IR infinite extra dimensions. 
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